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ABSTRACT: In this letter we propose to use an extension of the variational 
approach known as Truncated Conformal Space to compute numerically the 
Vacuum Expectation Values of the operators of a conformal held theory 
perturbed by a relevant operator. As an example we estimate the VEV’s of 
all (UV regular) primary operators of the Ising model and of some of the 
Tricritical Ising Model when perturbed by any choice of relevant primary 
operators. We compare onr results with some other independent predictions. 
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1 Introduction: the need of a variational approach for 

VEVs 


One physical application of (D = 2) conformal field theories, Jl[] , is that 
they can be successfully used to describe the critical properties of statisti¬ 
cal systems, such as critical indices and critical correlators, (see e.g. 0). 
Furthermore, following this line of thought, the behavior of the system near 
criticality, can be obtained by analyzing a conformal held theory perturbed 
by one or more relevant operators i.e. when to the conformal action Soft is 
added a perturbation of the form 

AS= Jdx'EWiM, (1.1) 


where A^ (<& is ) are bare couplings (operators) and 0 < Xi < 2 (aq = xo 
being the scale dimension of operator O). 


A wide class of perturbations have been found to generate integrable 
models, [j3|], for which the S matrix is known exactly and for which the off 
shell behavior of the theory can be estimated by use of the form factor method 
0, that gives rise to a long distance expansion for correlators. 


Another approach to the problem [|,i is obtained by using perturba¬ 
tive expansions around the (massless) conformal held theory, from which one 
could get informations on the short distance behavior of the complete the¬ 
ory. This point of view is important not only because it is complementary 
to the previous one, but also because it is in principle independent on the 
integrability of the model. The main problem in this case is the presence 
of I.R. divergences that arise in the naive Gell Mann Low series when the 
perturbing operator is relevant. These 1R problems can be solved by use 
of the Operator Product Expansion (OPE). The idea of this approach came 
out very early in the general context of quantum held theories and statistical 
mechanics [[5]] (see || for similar ideas in QCD), was hrst introduced in the 
context of perturbed conformal theories by Al. Zamolodchikov in || (see also 


TOj for an analog independent proposal), and was developed at all orders in 


a very general context in [|T2] . 


The idea of of the OPE approach is the fact that the so-called Wilson 
coefficients C^ b (r, A) that enter in the Operator Product Expansion for the 
complete theory, 


<*„(r)* s (0)>* ~ E^(O){*.(0))» (1-2) 

C 

are essentially short distance objects that can be taken to have a regular, 
I.R. safe, perturbative expansion with respect to the renormalized couplings 
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A*, while non analytic contributions to correlator can be confined in the (non 
perturbative) Vacuum Expectations Values (VEV’s). Main result of ]12| is 
an I.R. finite representation for the n th derivative of C£ b with respect to the 
couplings evaluated at A* = 0, involving integrals of (eventually renormal¬ 
ized) conformal correlators. The first order formulas have been applied to 
two models very interesting from the point of statistical mechanics and in¬ 
tegrate systems: the Ising model (IM) in external magnetic field 
the Tricritical Ising Model (TIM) at T > T c , [[14 
for earlier applications.) The results obtained were found to be in agreement 


13 


and 

(See also (50| and Jll| 


with form factor estimates [TR [IB], [T7[ in the intermediate region. 


It should be evident to the reader that the perturbative knowledge of 
Wilson coefficients is not sufficient to fully reconstruct correlators: actually 
one has to give some informations on VEV’s. 

The first, main information on VEV’s is obtained from the Renormaliza¬ 
tion Group equations for the operators to which we are interested: as well 
known in perturbative quantum field theories and from direct RG considera- 
TTJ , renormalized [] composite operators VEV’s evolve by RG 


tions, see e.g. 


as 


|([4>J> = r‘([J> 4 ]) 


(1.3) 


in which l is the logarithm of the scale and the matrix Y b a contains only pos¬ 
itive integer powers of the renormalized couplings A* and forbids the mixing 
with higher dimensional operators (it is lower diagonal if operators are or¬ 


dered with increasing dimensions). It is clear from Eq.( |1.3|) that when one 
coupling is considered, (as is the case for this paper) the mixing of two oper¬ 
ators, responsible of the presence of logarithms of A, could happen only if the 
difference of two operators is equal to a positive multiple of the dimension 
of A (resonance condition). When this is not the case the only nonperturba- 
tive information on operators VEV’s is contained in a constant, noted by A 
below: , 

([$ a ,])=A1(A',A)|A|^ (1.4) 

(A being the scale dimension of the perturbing operator). In the following, 
to avoid UV problems we will limit ourselves to compute the VEV’s of op¬ 
erators that do not need any multiplicative renormalization^ and for which 
the resonance condition does not hold. 


If the complete theory is integrable, the knowledge of the S matrix gives 
some informations on the bulk energy that can be recovered from the ex¬ 
pression of the free energy in terms of the mass scale as obtained from the 

1 We denote by [4>] any renormalized operator 4> 

2 This condition is satisfied for an operator <b A ' if A + A' < 1, as easily seen from an 
(IR regulated) perturbative expansion analysis for the VEV in question. 
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Thermodynamic Bethe Ansatz approach, |2l], . In this case (if the relation 

between the main mass and A is known, p3| , |22[| ) the VEV of the perturb¬ 


ing operator is simply the derivative with respect to A of the free energy. 
These informations are available for all the _Ad(£) + 01,3 models as well as all 
the _Ad(£) + 0i 2 (-M(£) + 02 i).[] In particular they were used in the papers 

BBSS' 

Moreover very recently Lukyanov and A. Zamolodchikov[|4|] gave an exact 
prediction for the VEV’s of the type (e*“^) in the Sine-Gordon model that, 
by quantum group reduction, was translated in the knowledge of all VEV’s 
of primary operators of all Ad(£) + 0i, 3 theories. 

In this letter we propose to extend the use of a variational approach, 
known as Truncated Conformal Space approach, first introduced by Al. 
Zamolodchikov in the context of perturbed conformal field theories a , to 
compute numerically the constants A above. The interest in our exten¬ 
sion is mainly that it completes the framework of short distance behavior 
of perturbed conformal field theories without requiring the integrability of 
the model. Moreover the estimates for TIM and IM that we give here as an 
example, constitute an independent check of previous evaluation of constants 
A and in particular an independent check of the Lukyanov Zamolodchikov 
expressions. 


In next section we review the variational method and the technical obser¬ 
vations of the present case. In section |3| we describe the numerical analysis 
and the results for the two examples comparing with known results, while 
conclusions are left to section [4|. 


2 The method 

To deal with the nonperturbative IR information, we need an IR regu¬ 
larization of the complete theory on the whole plane: to do that we will use 
the corresponding field theory on an infinite cylinder of circumference R and 
we will recover the original theory by taking the limit R —> oo. The regu¬ 
larized theory will be seen as a conformal theory perturbed by an operator 
$A (the tilde will refer to quantities of the regularized theory). The present 
approach assumes implicitly that the operators of the complete theories on 
the cylinder and on the plane are labeled by the operators of the unperturbed 
conformal theory, and that in the limit R —> oo states and operators on the 

3 We use the notation AtR), &i : k = d’A; k to denote respectively the minimal conformal 
theory with central charge c = 1 — and the scalar primary operator of scale dimension 

A i t k = and hs extension on the complete theory. 
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cylinder tend to corresponding quantities of the plane. The hamiltonian of 
the regularised theory will have the form: 


H = H 0 + V (2.1) 

where Hq is the Hamiltonian of the conformal theory on the cylinder and 
V = dw 2 &a(w) is the perturbation (w = W\ + iw 2 ). By use of the con- 

27T 

formal transformation z = e~R w the conformal theory on the cylinder can be 
shown to be equivalent to a theory on a punctured plane where z lives (this 
equivalence explains why in the limit R —> oo we recover the theory on the 
plane). In terms of quantities of the auxiliary punctured plane (denoted by 
an hat) we can write 


27r ~ - 

Hq — + Lq — c/12) 

(2,2) 

♦aW = \ 2 Z Z \ 2 a $a(z) 

Jr t 

(2.3) 


(we will fix without any loss of generality the quantization ”radial time” 
W\ = 0, corresponding to \z\ = 1). 

The key property that we will use is the fact that VEV’s of the complete 
theory on the plane can be recovered as a limit of the corresponding object 
on the cylinder: 

(0|$ A '|0) = lim (6|$ A '|0). (2.4) 

R —>oo 


Next step is to use the Truncated Conformal Space approach, to deal 
numerically with the theory at finite R, as first proposed in |20|. The TCS 
variational approach consists in truncating the conformal space of states 
(supposed to be the same for the complete theory) to some space of finite 
dimension N (obtained by giving an upper limit on values of x = L 0 + L 0 ) 
and to diagonalize the hamiltonian that will be described by a self adjoint 
matrix in this case. 

More in detail, denoting the set of (a priori nonorthogonal) conformal 
states with one must diagonalize the N x N hamiltonian: 


H) = g ik (k\H\j) (2.5) 

(i g lk above is the inverse of the matrix = (i\k) and is added to keep into 
account the nonorthogonality of the chosen basis that manifests also in the 
modified resolution of the identity 1 = \i)g l] {j\)- Remembering Eqs.( [2.1| - [2.3| ) 
it follows: 


rri _ 27T 

H >~n 


(Xi - c/\2)Sij + 2ttA(4L)2( i 

27T 


( 2 . 6 ) 
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The matrix elements (/c|<I>a|j) can be reconstructed from the knowledge of 
Wilson coefficients of the conformal theory on the punctured plane, see [|20| . 


TCS method has been used successfully to recover information on the 
spectrum of the complete theory on the cylinder by use of the secular equa¬ 


tions for Hy. 


det 


H) 


E^Si 


= 0. 


(2.7) 


It is clear that this method should converge for fixed R when N is big enough. 
When this condition is no more satisfied unphysical truncations errors are 
important and results are no more reliable. As suggested in |2(J a check 
on the behavior ~ R (when this is the case) in the R —> oo region by 
numerically verifying if 

xibiir - 1 (28) 

is a good control on the choices of R and N. 

The step forward that we propose here is to compute also eigenvectors 
il>W of Hj, in particular the ground state vector : from this information 
the VEV’s of the complete theory can be obtained by taking numerically the 
double limit R —> oo, N 


oo: 


($A')/I^| 1_A = lim (0|d > A'|0)Ar/|A|i-A 
R,N —>oc 


A(A',A) = 

(0|$A'|0)jV = 

(0|4>a'|o) k = (v , o ,v) ’( i l®^'b) A Ar)J ) / (i’o' 1 'aij’Po'y 


:^) 2A '(6|i>A-|o) K 


(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 


Another useful test on R , N (more sensitive to the VEV one is considering) 
is to verify if (0 | < Ia' |0)jv ~ R 2A ' (as imposed by the fact that A(A', A) is a 
finite number), i.e. to check numerically if 


1 dlog(0|l>A'|0)jv 

(2A') d log R 


= 1. 


( 2 . 12 ) 


3 Examples: IM and TIM 

In this section we will apply the ideas dclined in the previous section to 
the case of the Ising Model and Tricritical Ising Model (respectively Ad (3) 
and Ad(4) conformal theories). The primary scalar operators of the theory 
are (in the notation ($a, A)): 

IM = {(1,0), (a, 1/16), (e,l/2)} 
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±A, vac. 

1/16 

1/2 

+1/16, vac 

-1.277(2) 

1.94(6) 

+1/2, sym 

0 

U.V. 

—1/2, ssb 

±1.69(2) 

U.V. 


Tabic 1: Coefficients A a, a 1 for perturbed Ising Model (see Eq. (0)): A' is 
reported at the top of each column; the sign before A refers to the sign of 
the coupling; see text for explanations about the different vacua. 


TIM = {(1, 0), (a, 3/80), {£, 1/10), K, 7/16), ( £ ', 3/5), (/, 3/2)} 


The perturbation we will consider are cr, e for IM and cr, e, a', s' for TIM. 


For the numerical elaboration we used a slight modification of the Math- 
ematica program STRIP [19] that can deal minimal conformal field theories 
up to level 5 (i.e. up to N = 43 states for IM and N = 228 for TIM). The 
result of the analysis are reported in Table (1) for IM and in Table (2) for 
TIM. According to the type of the perturbations the spin reversal symmetry 
of the conformal theory can be exact, spontaneously broken or explicitly bro¬ 
ken in the complete model and we denoted the corresponding vacuum state 
(of the theory on the cylinder) respectively by ”sym, ssb, vac”. A detailed 
analysis of the structure of vacua in TIM can be found in |L8|. The values 
0 in tables are fixed by spin reversal symmetry, while with ”U.V.” we mean 
that the resonance condition is satisfied and that the VEV is no more of the 


form (|1.4p . 


The numerical strategy we used was (for a fixed value of N ) to look for 
the values of R such that Eq. (|2 .1 2|) is satisfied and to have a look at what 
happens at Eq.( [2.8p . The (indicative) errors reported are fixed by difference 
with the result at N = 8 for IS, N = 22 for TIM. The optimal situation (that 
we found in many cases) is that when Eq. ( |2 .1 2| ) is satisfied Eq. ( |2.8| ) holds 
within 1 — .1%. 


The general tendency we observe is that truncation errors are bigger for a 
perturbation <1 >a of higher A, and for fixed perturbation are bigger for VEV’s 
(<h A ') of higher A . 

More difficult (due to truncation errors) are the cases TIM + a' and 
TIM — s' 0 for which the truncation errors are bigger and stability is not 
achieved at the values of N considered. 

4 In TIM + s' all the VEV’s of relevant operators are 0 by symmetry apart from (e') 
that contains log’s. 
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±A, vac. 

3/80 

1/10 

7/16 

3/5 

+-^,vac 

-1.5396(8) 

1.336(3) 

-1.55(3) 

1.92(6) 

+Tn,sym 

0 

-1.466(7) 

0 

3.4(2) 

~Tn,ssb 

±1.594(2) 

1.466(4) 

±2.38(6) 

3.5(2) 


Tabic 2: Coefficients A a, a 1 for perturbed Tricritical Ising Model (see Eq.(|5])): 
A' is reported at the top of each column; the sign before A refers to the sign 
of the coupling; see text for explanations about the different vacua. 


In TIM + a' caseQ we can estimate (unsymmetric vacuum): 

.4(3/80, 7/16, +, 4) ~ -1.1 (3.1) 

4(1/10,7/16,+, 4) ~ 1.2 (3.2) 

while (e 1 ) contain log’s, and (a 1 ) is very unstable (this it not a surprise since 


it was observed 


that the level of accuracy of the STRIP program is 


not enough to describe the ground energy of which (a') is essentially the 
derivative) 

In the TIM — s' case we obtained 


±1.4 


(3.3) 


31(3/80,3/5,4) 
while (e) is not yet stabilized. 

We report now some predictions on the VEV’s that can be found in the 
literature (see footnote [5] for convention used). 

First of all we quote the exact predictions on the VEV of the perturbation 


coming from the Thermodynamic Bethe Ansatz ||2l|. [22]. [23 


4(1/16,1/16,+, 3) = -1.27758227605119295- 
4(1/10,1/10,+,4) = -1.468395424027621489- 


(3.4) 

(3.5) 


Secondly from the general formulas of Lukyanov and Zamolodchikov E4 
we can derive: 


4(1/16,1/2,-, 3) = 1.7085219053746979652-•• 
4(3/80, 3/5,-, 4) = 1.4105216239133927375 


(3.6) 

(3.7) 


We will use in the following the notations 


A / 


($ A d=4(A / ,A,SignA±)|A| — 


for the theory M. (£) + A<1>a ■ 



































We report then some results obtained by the authors in the previous 
analysis of critical Ising Model in magnetic held (from fit of lattice data and 
the imposition of C and A exact sum rules, [p^, p6]] ) 


.4(1/2,1/16, +, 3) = 2.02(10) 


(3.8) 


(we reported an error of 5% as suggested in the paper) as well as in TIM 
plus e (C and A exact sum rules imposed) 


4(3/5,1/10,+, 4) = 3.8(4). 


(3.9) 


The estimate 


4(1/2,1/16,+, 3) ~ 1.95 


(3.10) 


has been obtained recently by use of MonteCarlo simulations, [27 


ft appears that all the predictions are consistent with our results within 
the given errors! 


4 Conclusions 


The purpose of this letter was to suggest that a variational method could 
be efficient to recover (numerical) non perturbative informations on VEV’s 
independently of integrability of the underlying perturbed conformal theory. 
We want to emphasize in particular that such a method is in principle capable 
of giving informations on all VEV’s of the UV regular scalar operators of the 
complete theory, secondary operators and multiple couplings perturbations 
included. 


As an example of the application of the method we considered the per¬ 
turbed 1M and TIM, and we studied the VEV’s of primary operators of 
those models with a slight modification of an existing program [|lij, obtain¬ 
ing encouraging and nontrivial results already at level 5. Present results are 
in agreement within errors with exact results from TBA [pi], [23j] and from 


Lukyanov and A. Zamolodchikov [p[|], with some previous existing estimates 
of the authors, jT3], |T|] and with MonteCarlo numerical analysis |27J . Morever 


new original estimates have been given for some VEV’s, see tables. 

High precision estimates of constants 4, parameterizing UV regular VEV’s, 
are realizable with a more sophisticated (but not un-feasible) program to deal 
with conformal states and matrix elements, as explained in 
are outside the goal of this paper. 


19 and 
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We conclude by observing that the results of this letter should be consid¬ 
ered as an important step forward towards the practical description of short 
distance behavior of statistical mechanical systems in two dimension by use 
of the OPE approach: the all order formulas for the Wilson coefficients de¬ 
veloped in |L2| together with the present approach constitute a complete, 
consistent tool that can reach any accuracy, at least in the case of absence 
of UV divergences. 
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